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ABSTRACT 

Quasar-driven winds are currently the best candidates for accounting for the pre- 
heating of the intergalactic medium in clusters. Such winds, occurring during early 
phases of the evolution of spheroidal galaxies, shock-heat the interstellar gas, thus 
inducing a detectable Sunyaev-Zeldovich effect. We estimate the amplitude and the 
angular scale of such effect as well as its counts as a function of the comptoniza- 
tion parameter y. The contamination due to radio emission by the quasar itself is 
also discussed. The corresponding mean Compton distortion of the cosmic microwave 
background spectrum is found to be well below the COBE/FIRAS upper limit. 
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1 INTRODUCTION 



It is natural to expect that extremely powerful sources, such 
as quasars, strongly affect the surrounding medium. The re- 
cent evidences of a tight relationship between black hole 
masses and velocity dispersions of the host galactic bulges 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000a,b; Merritt 
& Ferrarese 2001a,b; McLure & Dunlop 2001) have strongly 
. ' strengthened the case for a close connection between the evo- 
, lutionary pathways of spheroidal galaxies and of quasars. 
5^ ■ Super-massive black holes have been found to be ubiqui- 
ty ' tous at the centers of local spheroidal galaxies (Kormendy & 
Richstone 1995; Magorrian et al. 1998; van der Marel 1999). 
The observed correlation between the mass of the black hole 
and that of the galaxy spheroidal component hints at a sub- 
stantial feedback of the nuclear activity to the surrounding 
medium (Silk & Rees 1998; Monaco et al. 2000; Granato et 
al. 2001): energetic quasar-driven winds can sweep up the 
interstellar gas and halt both star formation and the growth 
of the central black hole. 

If a significant fraction of the huge amount of energy 
released by quasars goes into ionization and heating of the 
neighboring gas, Compton cooling may produce a detectable 
Sunyaev-Zeldovich (1972; SZ) effect. The idea of strong 
shock heating of the medium originated by energetic out- 
flows from early quasars, originally developed by Ikeuchi 
(1981), was recently revived by Natarajan, Sigurdsson & Silk 
(1998), Natarajan & Sigurdsson (1999), and Aghanim, Bal- 
land & Silk (2000). Evidences of strong quasar-driven winds 
can be seen in Broad Absorption Lines (BAL) quasars, com- 
prising 10%-15% of optically selected quasars (Hamann & 
Ferland 1999 and reference therein). The dynamic interac- 



tion of such energetic outflows with the surrounding proto- 
galactic gas can heat it to high temperature. A release of 
an amount of mechanical energy from Active Galactic Nu- 
clei several times larger than that produced by supernovae in 
the surrounding galaxies may be necessary to account for the 
pre-heating of the intergalactic medium in clusters (Wu et al. 
2000; Kravtsov & Yepes 2000; Balogh et al. 2001; Bower et 
al. 2001). Natarajan & Sigurdsson (1999) suggested that the 
SZ effect associated to very energetic, quasar driven, winds 
may account for the reported isolated Cosmic Microwave 
Background (CMB) temperature decrements in directions 
where no clusters of galaxies are detected, but quasar pairs 
are present (Jones et al. 1997; Richards et al. 1997). De- 
tection of this effect would obviously be informative on the 
physics of quasar/galaxy evolution. 

In this paper we present a new investigation of the prob- 
lem, taking explicitly into account the relevant energetics as 
well as the observed relationships between black hole and 
host galaxy properties. Our approach, focussing on ener- 
getics, is less liable to uncertainties ensuing from the poor 
knowledge of the details of the gas heating process. In Sec- 
tion 2 we introduce the basic ingredients entering our calcu- 
lations and estimate the amplitude and the angular scale of 
the SZ effect. In Section 3 we present tentative estimates of 
its number counts. The results are discussed in Section 4. 



2 SZ EFFECT FROM QUASARS 

The main factors determining the amplitude of the quasar- 
driven SZ effect, usually measured in terms of the comp- 
tonization parameter 
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r R kTJr) 
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(where T e , n e and m e are the electron temperature, density 
and mass, cjt is the Thomson cross section, and R is the 
radius of the ionized region) can be illustrated as follows. 

The fractional amount of energy, Ae/ecMB, transferred 
to the CMB by Compton cooling of the plasma is related to 
y by (Zeldovich & Sunyaev 1969): 

1 Ae 
V--7- 



(2) 



4 ecMB 

where e CM B = aT^ MB ~ 4.2 x 10 _13 (1 + zf erg cm" 3 . 

In the redshift range where quasars are observed 
(z ,$ 6.3, Fan et al. 2001b) and for realistic values of the 
density of the medium, energetic cosmological blast waves 
remain adiabatic to relatively low z, i.e. the cooling time 
is longer than the expansion timescale, t cxp (Voit 1996; see 
also the discussion below). Thus, if a fraction fh of the total 
energy released by the quasar, E tot , goes into heating of the 
gas, the amount of energy per unit volume transferred to 
the CMB through comptonization in a time t < t c is: 

fhEtot t 



Ae: 



V t c 



(3) 



where V is the volume occupied by the hot gas and t c = 
3m e c/4o- T aTcMB - 7.3 x 10 19 (1 + z)~ 4 s is the Compton 
cooling timescale. We adopt, for simplicity, an Einstein- 
de Sitter cosmology, so that the expansion timescale is 
t cxp = a(t)/a(t) = (1/H )(l + z)~ 3/2 , where a(t) is the 
cosmic scale factor and 1/Ho — 6.17 x 10 17 h^g s, hso = 
iio/SOkms^Mpc- 1 . 

2.1 Evolution of the shock within the quasar host 
galaxy 

Close relationships between the mass of the central black 
hole, Mbh, and the properties of host elliptical galaxies or 
of galaxy bulges have been reported. Kormendy & Rich- 
stone (1995), Magorrian et al. (1998), and McLure & Dunlop 
(2001) found that Mbh is proportional to the luminosity of 
the spheroidal galaxy component. A tighter correlation be- 
tween Mbh and the line-of-sight velocity dispersion, a, was 
discovered by Gebhardt et al. (2000) and Ferrarese & Mer- 
ritt (2000), although the two groups find somewhat different 
slopes. We adopt the relationship obtained by Gebhardt et 
al. (2000), based on a larger sample: 

15/4 



Mbh = 1.2 10 8 



200 km/s 



M,-: 



(4) 



The total amount of energy released by the quasar is Etot = 
fflMbhC 2 , where e^ ~ 0.1 is the mass-to-energy conversion 
efficiency of the black hole. 

We model the host galaxy as an isothermal sphere with 
density profile 
-2 



p(r) -- 
with 



po 



(•») 



(6) 
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and an outer cut-off radius R g defined by the condition that 
the mean density within R g is S ~ 200 times higher than 



the mean density of the universe at the redshift Zf when the 
galaxy formed [numerical simulations (Cole & Lacey 1996) 
show that this radius approximately separates the virialized 
and infall regions; see also Navarro et al. 1997]: 



R 



■2a 



tfoS 1/2 (l + z/) 3/2 



V200 7 



(if:) 



kpc. (7) 



Several lines of evidence indicate that massive ellipticals 
were already in place by z ~ 2.5 (Renzini & Cimatti 1999; 
Ferguson et al. 2000; Daddi et al. 2000; Cohen 2001), al- 
though the issue is still somewhat controversial. For Zf ~ 2.5 
Eq. (fa) yields for R g values consistent with current estimates 
for massive galaxies (see, e.g., Peebles 1993). We assume that 
the distribution of baryons reflects that of dark matter. 

It is interesting to note, in passing, that this model, to- 
gether with the Gebhardt et al. (2000) relationship [Eq. (§)] 
implies Afbh oc M^ 5 , Mdm being the total mass of the 
galaxy, dominated by the dark matter halo. McLure & 
Dunlop (2001) find a best-fit relationship between black- 
hole mass and _R-band luminosity of the host galaxy bulge 
Mbh oc L^ 525 . Adopting, as they do, the relationship 
Mbuigo oc L 1 ' 31 , determined by j0rgensen et al. (1996), we 
have Mbh oc M^^ e . If we adopt the Magorrian et al. (1998) 
dependence, Mb u i gc oc L 1 , which according to Laor (2001) 
may be more appropriate, we get Mbh oc Mb u 2 g C . In both 
cases the observationally derived relationship between Mbh 
and Mbuigo is very close to the one yielded by the present 
model if Afbui gc oc Mum- 

The radius R s of a self-similar blast-wave carrying an 
energy Et = fhEtot, expanding in a medium whose density 
varies as r -2 , increases with time as (Ostriker & McKee 
1988): 



Rs = 



£Eb 
3po 



1/3 



,2/3 



(8) 



where £ w 1.5. 

Using Eqs. (Q) and (j|) it is easily checked that, for real- 
istic values of the parameters, the time, t g , required for the 
shock front to reach the outer radius of the galaxy 



8.9 x 10 16 h 



-3/2 / -Etot \ 1/ ' 6 / £fl 



(&)' 



Vio 62 / V0.1 
V200J 



-2/3 



-9/4 



(9) 



is always shorter than the expansion timescale t C x P . There- 
fore the shocks propagate outside of the host galaxy, and 
can heat up the general intra-cluster medium (ICM). How- 
ever, due to the much lower electron density in the ICM, 
compared to the mean density within the galaxy, most of 
the SZ signal comes from within R g (see also da Silva et al. 
2001). 

The main radiative cooling process for the redshifts of 
interest here (z < 6) is free-free, whose cooling timescale is: 



t« = 



3?7<g k'-L g 
wg 



(10) 



sy-j 1/2 2 1 c 

where w% = 1.410 T e n e g erg s is the cooling rate of 
a plasma with electron number density n e and temperature 
T e . The present mean ratio between baryon (mostly in stars) 
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and dark matter mass in massive spheroidal galaxies is esti- 
mated to be Mb/MnM ^ 0.03 (McKay et al. 2001; Marinoni 
& Hudson 2002). An upper limit to the cooling rate (and, 
correspondingly, a lower limit to tg) is obtained assuming 
that essentially all such baryons were in the interstellar gas, 
so that 



n e (r) = 3 x 10" 



2-irGm p r 2 



(11) 



where m p is the proton mass, we find (setting g — 1) that 
t B < t g [Eq. (§)] for 



r < r« 



«>,-W T - v 1/4 f E ^Y /20 f eR \ 
^ 50 Vio«kJ V 106 2 J Vol) 

■(&)' 



3/5 



1200/ 



(l + Z )- 9/8 kpc, (12) 



i.e. for only a very small fraction of the volume heated by 
the blast wave [see Eq. (^)]. Therefore we neglect, in the 
following, the radiative losses. 



2.2 Energetics of the shock 

The fraction of power released by the quasars which goes 
into heating of the gas is highly uncertain. Natarajan & Sig- 
urdsson (1999) assume this fraction to amount to about half 
of the bolometric luminosity. Analyses of the X-ray proper- 
ties of the intra-cluster medium (ICM) and, in particular, 
of the temperature-luminosity relation (Ponman et al. 1999; 
Valageas & Silk 1999; Cavaliere et al. 2000; Tozzi & Nor- 
man 2001; Balogh et al. 2001) indicate the need of a sub- 
stantial "pre-heating" of the gas before it collapses into the 
cluster. As discussed by Cavaliere & Menci (2001), mea- 
surements of the associated SZ effect very effectively probe 
the processes responsible for non-gravitational heating of 
the ICM. Bower et al. (2001) have shown that, unless the 
transfer of supernova energy to the ICM is very (perhaps 
unrealistically) efficient, additional heating sources, such as 
mechanical energy from quasar- driven winds are required. 
Adopting a baryon mass corresponding to fif, = 0.025/i~ 3//2 
(h — Ho/ 100) and a fraction of baryons converted into 
stars of /gai = 0.16/i 1//2 , Bower et al. (2001) estimate that 



esN — 1.3-2/i _1//2 = 1.8— 2.8h^Q^ 2 must go into heating of 
the ICM, and that most of it is unlikely to come from su- 
pernova explosions. Other energy sources, the most obvious 
being quasar winds, must be advocated. 

Adopting a ratio of /* = 0.03 between the mass in stars 
and the gravitational mass of spheroidal galaxies we have: 



h ~ 4.7xl0- 2 ^^(^) 
J 0.03 2.8 \0.lJ 



-4/5 



(E tot \ 
' \ 10 62 J 



-1/5 



3/2/ s y 
1 \2ooy 



V 3.5 ) 



-3/2 



(13) 



indicating that the mechanical power is a minor fraction of 
the bolometric luminosity. 



2.3 Amplitude and angular scale of the SZ effect 

As argued in sub-section 2.1, most of the SZ signal is ex- 
pected to come from within the galaxy, so that we will set, 
in Eq. (§), V = V g = (4n/3)R 3 g and t = t g . On the other 




Figure 1. Differential counts per deg 2 and per unit interval 
of log(j/) (the two upper curves) or of Iog(j/ O t, s ) (the two lower 
curves) of SZ signals due to quasar-driven blast waves, as a 
function of the comptonization parameter y (or j/ b s )i f° r the 
reference values of the parameters. The two solid lines corre- 
spond to t q = 10 7 yr, while the two dashed curves correspond 
to t q = 10 8 yr. The two lower curves take into account the dilu- 
tion effect of a Gaussian observing beam with FWHM = 10'. 



hand, only a minor fraction f g of the energy carried by the 
blast-wave goes into heating of the gas within the galaxy 
(most of it is dissipated in the ICM). Inserting in Eq. (El) 
the above results, using for fh the expression of Eq. (|l3|), 
and taking into account that, in the Rayleigh- Jeans region 
(AT/T) RJ = -2y we end up with 



2.2 x 10" 



(E tot \ i/15 
U0 62 / V2007 



4fg_ ftSN\ 1/2 /e R \ 
0.1 V 2.8 J VO.l/ 

f* x 1/2 



-4/15 



The angular radius is 8s z 



0.03 

R g /d A 



3/2 



cIa(z) — 



2c 1 



(1 + z 



1/2 



Ho 



(1 



(14) 



(15) 



being the angular diameter distance. From Eq. (bj) we have: 



Osz =s 17" 



,,/£t^\ 4/15 /£«\- 4/15 /l+£\ 
\W 62 ) \0.l) \ 3.5 J 



4 / 15 {l + z\-V 2 d A (2.5) 
<1a{z) 



,(16) 



where c£a (2.5) = 1595/ftso Mpc. The amplitude and the an- 
gular diameter of the effect, 29 sz — 34", for the reference 
values of the parameters, are not far from the results re- 
ported by Richards et al. (1997): \6T/T\ ~ W' 4 over an area 
of 30" x 65", and by Jones et al. (1997): \5T/T\ ~ 1.4 x 10 -4 
in a beam of 100" x 175". Since powerful high-z quasars are 
highly clustered (Croom et al. 2001), the larger than ex- 
pected observed angular scale might be interpreted as due 
to the combination of blast waves from neighbor quasars. 
It should be stressed, however, that in neither case the ex- 
planation in terms of quasar-driven blast waves is clearly 
supported by observations. A more conservative explana- 
tion is that the observed signals come from the SZ effect in 
previously unknown clusters. 
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3 COUNTS OF SZ SIGNALS FROM QUASARS 

The epoch-dependent "luminosity function" of SZ signals, 
4>sz(y,z) can be roughly estimated from the luminosity 
function of quasars, given the quasar lifetime t q . We have 
adopted one of the analytical evolutionary models for the 
B-band luminosity function, 4>b{Lb, z)d\og Lb dz, proposed 
by Pei (1995): the two-power-law model with /150 = 1, 
qo = 0.5, and optical spectral index a — 0.5 (f(v) oc u~ a ). 
This model works well up to z ~ 4.5, but appears to under- 
predict by a large factor the surface density of higher redshift 
quasars, which, however, are so rare (Fan et al. 2001a) that 
their contributions to the counts of the SZ signal is small. 

The B-band luminosity is related to the total energy 
released by E to t = fcsl/stq, where ks is the bolometric cor- 
rection for which we adopt the value fcs = 6 (the median 
value for the sample of Elvis et al. 1994, corrected to account 
for the different definition of Lb used by Pei 1995). If t q is 
independent of both luminosity and redshift, and consider- 
ing only the SZ effect within the host galaxy, we have: 

<t>sz{y,z) = ^<f> B {L B ,zYf , (17) 

where <f>sz(y, z) is the "luminosity function" of the SZ signal 
per unit log y and z intervals, and the factor 15/4 comes from 
d log Ls/d logy = d log E to t /d log y [Eq. Jl4|)], The quasar 
lifetime is still very uncertain. Recent estimates (Salucci et 
al. 1999; Monaco et al. 2000; Martini & Weinberg 2001) 
suggest t q to be in the range 8 x 10 6 -10 8 yr, values ~ few x 
10 7 yr being favoured (Haehnelt et al. 1998; Pentericci et al. 
2002; Ciotti et al. 2001). 

Inserting Eq. ( Jl3| ) in the expression for t q [Eq.(^)] we 
have: 



7.7 x 10 15 h 



3/4 /£tot\ 4/15 / eg 



V 1.0"- /' (o.l) 



-4/15 



\ 2.8 / V0.03/ 



-1/2 



(—) 

\2ooy 



~ 1/2 (l + z\ 
V 3.5 J 



-3/2 



(18) 



Although the amplitude of the SZ effect is distance indepen- 
dent, the observed signal is affected by distance dependent 
beam dilution or resolution effects, so that the SZ signal 
observed with an instrumental solid angle Wbeam is ; 



y obs ~ yF(cj bc!l m/uJsz) , 
where iqsz = tt#sz and 



F(uj hea . m /ujsz) 



/losz if 

? ±8/lB 



(19) 



(20) 



For given z, F oc E tot oc y , the sign of the exponent 



depending on the ratio Wbeam/wsz [see 



Eq. ©] 



so that 



J/obs oc y 
then write: 

dNsz 



dy bs ln(10) y a b 



The differential counts per steradian of y a 



, dV , . , _, . dlogw ,„„. 
dz—<j>sz {yobs F,z) &y , (21) 
dz d log j/ b s 



where, for a Friedman universe, the volume element per unit 
solid angle and unit dz is dV/dz = (c/H )d 2 L (l + z)~ 13/2 , 
dL = (2c/ Ho)(l + z — y/1 + z) being the luminosity distance. 
The differential counts dNsz /dlogy = y\n(10)dNsz /dy are 
shown in Fig. 1 for two values of t q , 10 7 yr and 10 s yr. For 
both cases, the effect of dilution by a Gaussian beam with 
FWHM = 10' is also shown. 



The extremely steep slope of the bright end of the 
counts is model dependent. It comes out from the combined 
effect of the steep slope of the high luminosity portion of the 
quasar luminosity function and of the mild dependence of y 
on Etot [see Eq. (|l4|)] , implied by the present model. Clearly, 
the assumption of dispersion-less relationships between the 
quasar blue luminosity and Etot, and between Etot and y is 
an over-simplification. Dispersions will result in a flattening 
of the bright end of the counts. 

The surface density of SZ signals is strongly depen- 
dent on the amount of thermal energy injected in the in- 
terstellar gas of the host galaxy, which, according to the 
present model, is proportional to t^ 5 . For t q in the range 
10 7 -10 8 yr, we expect some 10 3 SZ signals per square de- 
gree with y > 2 x 10 -4 to be detected with an imaging 
survey with £ 30" resolution (FWHM). For higher values 

of y, the expected surface density is strongly dependent on 
the quasar lifetime t q . The survey should be carried out at 
10 GHz < v < 150 GHz to avoid blurring of the SZ signal by 
local radio or dust emission. The signal is strongly diluted in 
most current surveys aimed at mapping CMB anisotropies 
with angular resolution ^ 5'-10'. Even high sensitivity, all 
sky surveys such as those to be carried out by ESA's Planck 
mission, cannot efficiently detect such signals. Much better 
prospects for detecting SZ effects from quasar-driven blast 
waves are offered by ground based bolometric arrays and 
interferometers with arc-min or sub-arcmin resolution, such 
as the Arcminute MicroKelving Imager (Kneissl et al. 2001), 
and AMiBA (Lo et al. 2001). 

Since, in the present framework, quasars are signposts of 
massive galaxies at high redshifts and are therefore strongly 
biased tracers of the matter distribution, we expect them 
(and the associated SZ signals) to be highly clustered. 



4 DISCUSSION 

4.1 Effect of local radio emission 

It is entirely plausible, and perhaps required by evidences of 
a substantial preheating of the intra-cluster gas, that quasars 
inject in the surrounding medium an amount of energy suf- 
ficient to produce a detectable SZ effect. The corresponding 
angular scale is expected to be relatively small (sub-arcmin), 
so that the blurring effect by the emission from the quasar 
itself or from the host galaxy may be important. 

To estimate to what extent the radio emission associ- 
ated with the quasar might blur the SZ signal, we refer to the 
median ratio of monochromatic luminosities at 5 GHz and 
at vb — 6.82 x 10 14 Hz (corresponding to As = 0.44 fim), for 
the radio-loud and radio-quiet quasars in Table 2 of Elvis et 
al. (1994). We find log (ZsgW'b) median = -0.47 for radio- 
quiet quasars and = 2.2 for radio loud. 

Taking into account that the contribution to the an- 
tenna temperature at the frequency v within a solid angle 
(j->sz = tt^Iz, of a source of 5 GHz flux Ssghz and spectral 
index a ~ 0.7 (S v oc u~ a ) is: 



AT A 



S 5G Hz(v/5GR Z )- a c 2 



(22) 



2kbV 2 uJsz 

we have, for radio-quiet quasars, in the Rayleigh- Jeans re- 
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gion 

AT A 
T A 



3 (Ei2L\ V15 f?*y 



L11 °" v 10- y 

10 7 yr / 5GHz N^ ' 

to \ V ) 



(23) 



In the case of radio loud quasars the coefficient would be 
about 500 times higher and therefore the radio emission 
would easily overwhelm the SZ signal, except, perhaps, at 
high frequencies (~ 100 GHz). However, only a minor frac- 
tion ( £ 10%) of quasars are radio loud. But even in the case 
of radio quiet quasars the radio emission may, at least par- 
tially, fill up the SZ dip at v £ 10 GHz, particularly if the 

quasar is surrounded by an intense starburst, which is also 
radio bright, as is frequently the case at high z (Omont et 
al. 2001). The contamination by radio emission sinks down 
rapidly with increasing frequency, but at v Si 150 GHz, dust 
emission powered by star formation in the host galaxy, may 
take over, due to its spectrum steeply rising with increasing 
frequency. 

On the other hand, the quasar lifetime and the duration 
of intense starbursts are generally shorter than the time t g 
for the blast-wave to reach the boundary of the host galaxy. 
Therefore the SZ signals may be observed when the quasars 
that originated them are dead. 



4.2 Global comptonization distortion 

The mean distortion of the CMB spectrum due to comp- 
tonization by the electrons heated by the blast waves, mea- 
sured by the mean value (y) of the parameter y integrated 
along the line of sight, is 1/4 of the fractional amount of 
energy injected in the CMB per unit volume: 



iv) = 



Jh,eB 



dz d log L, 



<$>B {Lb , z)kBL/B £ exp dt 



£cmb(^) 



t c dz 



(24) 



where /h, c ff is the effective fraction of the total energy re- 
leased by the quasars carried by blast waves. From the above 
discussion it follows that pre-heating of the ICM requires 
/h.cff ~ 0.1. Using the double-power-law quasar luminosity 

function by Pei (1995), we find: 



(y) ~ 2.4 x 10 



-6 /h,. 



0.1 



(25) 



Using Pei's (1995) exponential luminosity function, the 
value of (y) increases by 20%. On the other hand, it may be 
noted that Eq. ( p4| ) includes the contribution (which turns 
out to be quite substantial) of low-z quasars to the integral 
(the contribution from the redshift range 0-2 is comparable 
to that from 2 < z < 4); however, in the present frame- 
work, the most energetic blast-waves should be associated 
to early phases of the quasar/galaxy evolution. Thus the ex- 
pected global comptonization distortion of the CMB spec- 
trum induced by quasar driven blast waves is well below the 
COBE/FIRAS limit (y) < 1.5 x 10" 5 (Fixsen et al. 1996). 
This limit could be (marginally) exceeded only in the unre- 
alistic case that the total energy emitted by all quasars over 
the entire life of the universe went into heating of the sur- 
rounding gas. This would also entail an excessive (by almost 
one order of magnitude) pre-heating of the ICM. 
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